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One loop mediated charged Higgs bosons decays — » W^V, V = Z, 7 are studied in 
the Minimal Supersymmetric Standard Model (MSSM) with and without CP violating 
phases. We evaluate the MSSM contributions to these processes taking into account 
B — t X B j constraint as well as experimental constraints on the MSSM parameters. In 
the MSSM, we found that in the intermediate range of tan [3 <J 10 and for large At and 
large fi, where the lightest top squark becomes very light and hence non-decoupled, the 
branching ratio of H — » W^Z can be of the order 10 — 3 while the branching ratio of 
H ± -+ W ± 7 is of the order 10 5 . We found also that the CP violating phases of soft 
SUSY parameters can modify the branching ratio by about one order of magnitude. We 
also show that MSSM with CP violating phases lead to CP-violating asymmetry in the 
decays — > W~*~V and H~ — » W~V. Such CP asymmetry can be rather large and 
can reach 80% in some region of parameter space. 



1. Introduction 

Supersymmetric (SUSY) theories, in particular the Minimal Supersymmetric Stan- 
dard Model (MSSM), are currently considered as the most theoretically well mo- 
tivated extensions of the Standard Model. Recently, phenomenology of the MSSM 
with complex SUSY parameters has received growing attention LU. Such complex 
phases provide new sources of CP violation which may explain electroweak baryo- 
genesis scenarios and CP violating phenomena in K and B decays It has 
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been shown in ^ that by assuming universality ol the gaugino masses at a high 
energy scale, the effects of complex soft SUSY parameters in the MSSM can be 
parameterized by two independent CP violating phases: the phase of the Higgsino 
mass term /i (Arg(/x)) and the phase of the trilinear scalar coupling parameters 
A = Af (Arg(Af)) of the sfermions /. The presence of large SUSY phases can give 
contributions to electric dipole moments of the electron and neutron (EDM) which 
exceed the experimental upper bounds. In a variety of SUSY models such phases 
turn out to be severely restricted by such constraints, i.e. Arg(^i) < (10~ 2 ) for a 
SUSY mass scale of the order of few hundred GeV ^. However, the possibility of 
having large CP violating phases can still be consistent with experimental data in 
any of the following three scenarios: i) Effective SUSY models , ii) Cancellation 
mechanism^ and iii) Non- universality of trilinear couplings A/^. In the MSSM, 
after electroweak symmetry breaking we are left with 5 physical Higgs particles 
(2 charged Higgs H^, 2 CP-even H°, h° and one CP-odd A ). In this study, our 
concerns is the charged Higgs decays H — > W V, V = Z, 7 with and without 
CP violating phases. Those channels have a very clear signature and might emerge 
easily at future colliders. For instance, if — » W^Z is enhanced enough, this 
decay may lead to three leptons final state if both W and Z decay leptonically and 
that would be the corresponding golden mode for charged Higgs boson. Charged 
Higgs phenomenology has been extensively studied in the literature. It has been 
shown that SUSY one-loop contributions can lead to decay rates asymmetry of 
H ± — > tb, bi, H ± — > Uid*,u*dj and H ± — > xfx^-iX^Xj ■ Similar study has been 
done for the single charged Higgs production at hadron collider ^ and shown that 
SUSY CP phases can lead to CP asymmetry. 

In this paper, we will discuss both the branching ratios of — > W^V with 
V = Z, 7^^as well as the CP-violating asymmetry in the decay rates H + — > W + V 
and H~ — > W~V that emerge from the presence of CP violating phases in MSSM. 

The paper is organized as follows. In section II, we describe our calculations and 
the one- loop renormalization scheme we will use for — > W^V and present the 
source of CP violation in H — * W^V. In Section III, we present our numerical 
results and discussions, and section IV contains our conclusions. 



2. CP violation in Charged Higgs decay: H^ 1 — ► VF^V 
2.1. One loop amplitude H^ 1 — > W^V 

In MSSM, at tree level, the coupling _ff ± VU ± 7 is absent because of electromagnetic 
gauge invariance t/(l) em , while the absence of H^W^Z is due to the isospin sym- 
metry of the kinetic Lagrangian of the Higgs doublet fields. We emphasize here 
that it is possible to construct models with an even larger scalar sector than 2 
Hig gs doublets. One of the most popular being the Higgs Triplet Model (HTM) 
A noteworthy difference between MSSM and HTM is that the HTM contains a 
tree level ZW^H^ coupling. Therefore, in MSSM, decays modes like if ± — > VF ± 7, 
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H ± -> W ± Z are mediated at one loop level I 12 | 13 | 14 | 15 | 16 | Although these de- 
cays are rare processes, loop or/and threshold effects can give rise to substantial 
enhancement. Moreover, once worked out, any experimental deviation from the re- 
sults within such a model should bring some fruitful information on the new physics 
and also allow to distinguish between models. Let us give the general structure of 
the one loop amplitude for — > W^V and discuss the renormalization scheme 
introduced to deal with tadpoles and vector boson-scalar boson mixing. 

The amplitude M. for a scalar decaying to two gauges bosons V\ and V2 can be 
written as 

a e,/ 

where ey ; are the polarization vectors of the Vi. According to Lorentz invariance, 
the general structure of the one loop amplitude of S — > V^V^ decay is 

(2) 

where pi ; 2 are the momentum of V\, V% vector bosons, ^1,2,3 are form factors, 
and e^pa is the totally antisymmetric tensor. The form factor T\ has dimen- 
sion 2 while the other are dimensionless. Therefore, it is expected that in case of 
H ± —* W^Z decay, T\ will grow quadratically with internal top quark mass while 
^2,3 w iU have only logarithmic dependence Quite contrary, for i? ± — * W ±r f 
decay, the electromagnetic gauge invariance implies that T\ and Ti are related 



by T\ =l/2(m| l/ — rnj I ±)3 : 2 This means that only T2 and Tz will contribute 
to the decay H ± — > W ±n f and then the amplitude of — > W ±r y will not grow 
quadratically with internal fermionic masses. In case of if ± — ► W^Z, there is no 
such constraint on the form factors. In fact the CP asymmetry in this process 
H ± — > W^Z is very large compared to the CP asymmetry in H + — > i& The 
reason is that iJ* — > W^Z is loop dominated, while i? + — > i6 is dominated by the 
tree level contribution which is CP conserving. 

The one- loop amplitude for — ► W^V is evaluated in the 't Hooft-Feynman 
gauge using dimensional regularization. The typical Feynman diagrams that con- 
tribute to if± -y W±V are depicted in Fig. 1. Those diagrams contain vertex 
diagrams (Fig. 1.1 -> 1.11), W ± -H ± mixing (Fig. 1.12 -> 1.14), J ff ± -G ± mixing 
(Fig. 1.15 -> 1.17) and H ± -W ± mixing (Fig. 1.18 -> 1.20). 

Note that the mixing J ff ± -VF ± (Fig .1.12, 1.13, 1.14) vanishes for an on-shell 
transverse W gauge boson. There is no contribution from the W ± -G^ mixing 
because jG ± H T and ZG ± H T vertices are absent at the tree level. All the Feynman 
diagrams have been generated and computed usin g FeynArts and FormCalc L^TJ 
packages. We also used the fortran FF-package in the numerical analysis. 

Although the amplitude for our process is absent at the tree level, complica- 
tions like tadpole contributions and vector boson-scalar boson mixing require a 



a As it has been shown in Ref. |13j . the top quark contribution does not decouple while squarks 
contributions does. 
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Fig. 1. Generic contributions to H — > W^V 



careful treatment of renormalization. We adopt, hereafter, the on-shell renormal- 
ization scheme of Ref. Q]5], for the Higgs sector. In this scheme, field renormalization 
is performed in the manifest-symmetric version of the Lagrangian. A field renor- 
malization constant is assigned to each Higgs doublet $1,2- We follow closely the 
same approach adopted in Refs. [10l [2Qj for W^-H^ and H ± -G T mixing as well 
as renormalization of - 11 : // ' ;ind 7,W ± H T . We will use the following on shell 
renormalization conditions I 

• The renormalized tadpoles, i.e. the sum of tadpole diagrams T^^h and 
tadpole counter-terms 8h,H vanish: 

T h + 5t h = 0, T H + 5t H = . 

These conditions guarantee that the vevs i»i,2 appearing in the renormalized 



1, 2008 18:36 WSPC/Guidelines-IJMPA cpt'f 



H 



CP violation in Charged Higgs Bosons decays H+ — » W+ (7, Z) in MSSM 5 

b,X bj s b,x 



t,x 

2.1a 



W 




H " 



\ W 



t, h 

2.1b 



\ w 



v 

2.1c 



H 




t,X 
2.1d 



G 



H " 



6,/ 

i, ^ 
2.1e 



\ G 



H 



V 

2.1f 



G,H 



H 




2.2a 



H 




2.2b 



► -X ► 

H G 




W 




H W 



V 



2.2c 



2.2d 



Fig. 2. Generic contributions to H ± -> W± and H± -> G± mixing as well as counter-terms 
needed. 



Lagrangian Cn are located at the minimum of the one-loop potential. 
• The real part of the renormalized non-diagonal self-energy S#± y/± (k 2 ) van- 
ishes for an on-shell charged Higgs boson. This renormalization condition 
determines the counter term for 5[A V W+H T ] and S[Z„W+H T ]. 

Using the Slavnov-Taylor identity, which relates H ± W ± and H ± G ± self energies, 
the last renormalization condition about H+W+ mixing is sufficient to discard the 
real part of the H -G ± mixing contribution as well. 



2.2. CP violation in H± —> W ± V 

The CP violating phases of MSSM can lead to CP-violating asymmetry in the decay 
rates H + — > W + V and H~ — > W~V. In what follow, we will study the following 
CP-violating asymmetry 



CP _ T(H+ -> W+V) - T(H~ -» VK-y) 

~ r(ff+ -> VF+y) + r(#- -> w-u) 



(3) 



and show that the MSSM with complex phases together with the absorptive part 
emerging from some loop integrals can lead to non- vanishing 6 CP . Gcnerically, in 
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the case of — > W 7, the form factors J-\^ take the following form 
^ 3 (F+ - W + 1 ) oc J2 blwFViM + bl ijk PV 2 , ijk 

i,j,k 

where (ft) &) a ,yfc ( a — 1)2 and z,j and k label the particles inside the loop) are 
combination of MSSM couplings, namely charged Higgs couplings to charginos- 
neutralinos and charged Higgs couplings to squarks. PV a ,ijk and PY a ,ijk are Passarino 
Veltman functions. The CP asymmetry takes the following form 

S CP oc (|J-i(i/-)| 2 - \T 3 (H-)\ 2 - \Fi{H+)\ 2 + \T 3 (H+)\ 2 )/ a ll 

K E E ^(^,«feOftijfc)9*n(PV„ lii fcPV^ iiifc ) + (a <— > /3,PV — PV) 

It is clear from the above equation that in order to have a non-vanishing 8 CP we 
need both complex couplings Cbot,ijk & s well an absorptive part from the Passarino- 
Veltman functions. For example, diagrams Fig. 1.1c and Fig. 1.18 (resp. Fig. 1.2 and 
Fig. 1.19) will develop some absorptive parts if the decay channels H ± — ► xfx'j 

(resp. — > q^q'j) is open, i.e. niH± > m x m x ( res P- m H± > m ^ m ~r .)■ 
3. Numerics and discussions 

In our numerical evaluations, we use the following experimental input quantities ^4 
or 1 = 129, m z , m w , m t , m b = 91.1875, 80.45, 174.3, 4.7 GeV. In the MSSM, we 
specify the free parameters that will be used as follow: i) The MSSM Higgs sector 
is parameterized by the CP-odd mass m^o and tan/3, taking into account one- 
loop radiative corrections from [22] , and we assume tan/3 > 3. ii) The chargino- 
neutralino sector can be parameterized by the gaugino-mass terms Mi, M%, and the 
Higgsino-mass term \i. For simplification, GUT relation M\ ss M 2 /2 is assumed, in) 
Sfermions are characterized by a common soft-breaking sfermion mass M$usy = 
Ml = Mr, (i the parameter and the soft trilinear couplings for third generation 
scalar fermions A ti b.r- For simplicity, we will take A t — Ab = A T . 

When varying the MSSM parameters, we take into account also the following 
constraints: i) The extra contributions to the dp parameter from the Higgs scalars 
should not exceed the current limits from precision measurements 1 ^* \op\ < 0.003. 
ii) b — > sj constraint. The present world average for inclusive b — > 57 rate is 
B(B -» X sl ) = (3.3 ± 0.4) x 10~ 4 . We keep the B X sl branching ratio in the 
3cr range of (2.1-4.5)xl0~ 4 . The SM part of B -> X sl is calculated up to NLO 
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Fig. 3. Branching ratios of — * W^Z (left) and — * W^f (right) as a function of m H ± 
in the MSSM and 2HDM for M SUSY = 500 GeV, M 2 = 175 GeV, n = -1.4 TeV and At = A b = 
A T = —/Li for various values of tan (3. 



using the expression given in 23J. While for the MSSM part, the Wilson coefficients 
C7 and Cg are included at LO in the framework of MSSM with CKM as the only 
source of flavor violation and are taken from [24]. Hi) We will assume that all SUSY 
particles sfermions and charginos are heavier than about 100 GeV; for the light CP 
even Higgs we assume m^o > 98 GeV and tan [3 > 3 ^51. For charged Higgs boson 
we assume that mn± > 78 GeV from LEP experiments ^21. As the experimental 
bound on m h a is concerned, care has to be taken. Since we are using only one- 
loop appro xim ation for the Higgs spectrum, and as is well known, higher order 
corrections may reduce the light CP-even Higgs mass in some cases. It may be 
possible that some parameter space region, shown in this analysis, which survives 
the experimental limit m^o > 98 GeV with one loop calculation may disappear, 
once the higher order correction to the Higgs spectrum are included. 

The total width of the charged Higgs is computed at tree level from [28] without 
any QCD improvement for its fermionic decays —>//'. The SUSY channels like 
H + — > fifj and H + — » XiX^ are included when kinematically allowed. In Fig. [3] we 
show branching ratios of — > W^Z (left) and H ± — ► W^j (right) as a function 
of charged Higgs mass for tan/3 = 16 and 25. In those plots, we have shown both 
the pure 2HDm[B and the full MSSM contributions. As it can be seen from those 
plots, both for F ± -> W ± Z and ff ± -» 1V ± 7 the 2HDM contribution is rather 
small. Once we include the SUSY particles, we can see that the branching fractions 
get enhanced and can reach 10 -3 in the case of i/ ± W ± Z and 10" 5 in the case 
of -> W ± j. The source of this enhancement is mainly due to the presence of 



b Pure 2HDM means that we include just the 2HDM part of the MSSM that contributes here in 
the loop, i.e. only SM fermions, gauge bosons and Higgs bosons with MSSM sum rules for the 
Higgs sector. 
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Fig. 4. Branching ratios for — > W^Z (left) and — > W "f (right) as a function of At in 
the MSSM with M SUSY = 500 GeV, M 2 = 200 GeV, m H ± = 500 GeV, At = A b = A T = 
and — 2TeV < fi < —0.1 TeV for various values of tan/3. 



scalar fermions contributions in the loop which are amplified by threshold effects 
from the opening of the decay — > Ubj. It turns out that the contribution of 
charginos neutralinos loops do not enhance the branching fractions significantly as 
compared to scalar fermions loops. The plots also show that, the branching fractions 
are more important for intermediate tan/3 = 16 and are slightly reduced for larger 
tan/3 = 25. 

We now illustrate in Fig. [4] the branching fraction of _ff ± — > W^Z (left) and 
H± _» VF ± 7 (right) as a function of A t = A b = A T = -/i for M S usy = 500 
GeV and M-i — 200 GeV. Since b — > sj favors A t and /i to have opposite sign, 
we fix /i = —At and in this sense fi is also varied when A t is varied. Both for 
#± w ± Z and i/ ± W ± 7, the char gino-ncutralino contributions, which are 
rather small, decrease with /i = — A t : the larger At, the smaller chargino-neutralino 
contributions. As one can see from these figures, the plots stops at A t — 1.1 TeV and 
tan/3 = 3, because for larger A t the dp constraint will be violated. For tan/3 = 10 
and 20, the plots stop for the same reason. In case of — > W ± Z, for A t < 1 TeV 
it is the pure 2HDM contributions which dominate and that is why it is almost 
independent of A t , while for larger A t the branching ratios increase with A t . It is 
clear that the larger A t is the larger the branching ratios, which can be of the order 
of 10 -3 for — > W^Z with tan/3 = 10, for example. As we know from h° — > 77 
and h° — > 7Z in MSSM the squarks contributions decouple except in the light 
stop mass and large A t limit ^QJ i n JJ^ — * W^V case, the same situation happens. 
As we can see from Fig. |4] (left), for intermediate A t , 300 < A t < 1000 GeV, 
the squarks are rather heavy and hence their contributions are small compared to 
2HDM one. While for large A t the stop becomes very light < 200 GeV and hence it 
can enhance the — > W^V width. Of course this enhancement is also amplified 
by H^tL^pbji L couplings which are directly proportional to A t: b- In the case of 



February 1, 2008 18:36 WSPC/Guidelines-IJMPA cpt'f 



CP violation in Charged Higgs Bosons decays H^ 1 —> (7, Z) in MSSM 9 




p 10- 

+ 

t 
+ 
X 

m i - f 



tan|3 




tan(5 = 


= 10 


tan(5 = 


= 20 







-20 20 
Arg(A,)(deg) 



-20 20 
Arg(A t )(deg) 



Fig. 5. Branching ratios for H± -* W ± Z (left) and -+ W /± 7 (right) in the MSSM as a 
function of Arg(A t ) : M SUSY = 500 GeV, M 2 = 150 GeV, m H ± = 500 GeV, A t = A b = A T = 
—li = l TeV and for various values of tan 0. 

— > M /± 7 decay, the pure 2HDM and sfermions contribution are of comparable 
size, the branching ratio increases with A t . We have also studied the effect of the 
MSSM CP violating phases on charged Higgs decays. For H ± — > W^Z and i? ± — > 
W ±n j decays which are sensitive to MSSM CP violating phases through squarks 
and charginos-neutralinos contributions, it turns out that the effect of MSSM CP 
violating phases is important and can enhance the rate by about one order of 
magnitude. For illustration we show in Fig. [5] the effect of A t ^,r CP violating 
phases for M SU sy = 500 GeV, A tAr = -fx = 1 TeV and M 2 = 150 GeV. To avoid 
conflict with EDM constraint we set the /i parameter to be real^. As it is clear, the 
CP phase of A t: b,r can enhance the rates of both i? ± — » W ± {Z 1 7) by more than 
one order of magnitude. The observed cuts in the plot are due to b — > sj constraint. 

In Fig.© we show the CP asymmetry Sep for H ± — > W Z (left) and J? — > 
W ± "f (right) as a function of charged Higgs mass. As one can see in both plots, 
for rriH± 5s 512 GeV, Sep is very small, of the order 10~ 3 . In this case, the main 
contributions to Sep arise both from diagrams Fig. 1.1 and Fig. 1.18 with absorptive 
parts and real couplings and diagrams Fig. 1.2 and Fig. 1.19 with complex couplings 
and no absorptive part (the channel H — > tib* is not yet open!). However, once 
mH± <^ 512 GeV, the channel i? ± — > tib* is open, the diagrams Fig. 1.2 and Fig. 1.19 
have an absorptive part and consequently they contribute to Sep, which can go up 
to several percent. The thresholds of H ± — ► t\b\ at m#± w 512 GeV and of 
H ± —* tj)*, at m H± « 740 GeV are clearly visible in Fig.©. 



4. Conclusion 

In the framework of MSSM, we have studied charged Higgs decays into a pair 
of gauge bosons, namely, H*- — > W^Z and i? ± — » W^^f. We have also studied 
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the effects of MSSM CP violating phases in these processes. In contrast to previous 
studies, we have performed the calculation in the 't Hooft-Feynman gauge and used 
a renormalization prescription to deal with tadpoles, W^-H^* 1 and G —H mixings. 
The study has been carried out taking into account the experimental constraints 
on the p parameter as well as the inclusive 6-decay b — > sj. Numerical results for 
the branching ratios have been presented. In the MSSM, we have shown that the 
branching ratio of iJ ± — > W^Z can reach 10~ 3 in some cases while — ► W^j 
never exceed 10 -5 . Branching ratio of the order 10 -3 might provide an opportunity 
to search for a charged Higgs boson at the LHC through H^- — > W Z. The effect of 
MSSM CP violating phases is also found to be important and can change the size of 
the branching ratios of H ± — > W^V by about an order of magnitude compared to 
the CP conserving case. The CP-violating asymmetry in the decays H + — > W + V 
and H~ — > W~V can be rather important and can reach 80% in some region of 
parameter space [31] . 
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